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Abstract 

We prove existence of an optimal transport map in the Monge- 
Kantorovich problem associated to a cost c{x, y) which is not finite 
everywhere, but coincides with \x — yp if the displacement y — x be- 
longs to a given convex set C and it is +00 otherwise. The result is 
proven for C satisfying some technical assumptions allowing any con- 
vex body in R-^ and any convex polyhedron in W^, d > 2. The tools 
are inspired by the recent Champion-DePascale-Juutinen technique. 
Their idea, based on density points and avoiding disintegrations and 
dual formulations, allowed to deal with L°° problems and, later on, 
with the Monge problem for arbitrary norms. 
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1 Introduction 

The optimal transport problem, introduced by G. Monge in [H] at the end 
of the 18th century, has become very famous in the past 20 years, when the 
relaxed formulation given by L.V. Kantorovich in the '40s, coupled with the 
most recent advances, has finally allowed to give a complete solution to many 
useful issues of the problem. Given /q « dx and fi two probability measures 
on M.'^, the issue is to push forward /o to /i with a map T : M.'^ ^ M.'^ called 
optimal transport map which solves the following optimization problem: 



where c : [0, +oo] is a fixed cost function. The original problem 

by G. Monge concerned the case where c is the Euclidian distance, but later 
on many other cost functions have been investigated, and in particular the 
quadratic one c{x,y) = \x — which has been popularized by [5J. We will 
not enter now into the details of the formulation that L.V. Kantorovich gave 
to overcome the difficulties of the nonlinear behavior of Problem (M) since 
we want to present as soon as possible the particular cost function that we 
want to consider. In this paper we will actually focus on a quadratic cost 
with a closed convex constraint C gM.'^: 



where we denote by | ■ | the Euclidian norm. 

Solving (M) for this cost seems a natural question, both for its modelization 
meaning (imagine that C is the unit ball in W^: in such a case the problem 
aims at finding the optimal displacement minimizing the total effort, subject 
to a pointwise constraint \T{x) — x| < 1) and for mathematical applications. 
Actually, it appears in [7] that a strategy to solve more general transport 
problems passes through the solutions of some constrained cases, where it 
is imposed that T(x) — x belongs to a given set, which is often the face 
of a convex body. The L°° problems studied in [10] can also be considered 
constrained problems (if the minimal L°° norm | \T{x) — X 1 1 i^oo is L, this means 
that we could consider maps constrained to satisfy T{x) — x E B{0,L)). 
In many of these cases, once the problem is transformed into a convexly 
constrained one, then it is only necessary to find an arbitrary transport map 
satisfying such a constraint. The minimization of the quadratic cost among 
them is just an additional criterion to select a special one, and it is chosen 
since it usually guarantees better properties. We refer anyway to [7] for a 



(M) 





y — if y — X E C 
■foo otherwise. 



(1.1) 
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more detailed discussion about convex costs in general and quadratic costs 
with convex constraints. 

To illustrate the difference between our setting and the usual quadratic 
case c{x,y) = \x — yp, we give an example. 

Exemple 1.1. Let Ui, with i = 1, 2, 3 three points of the plan with the following 
distances: \ui — U2\ = \u2 — M3I = 1 and \ui — u^l = 5/4. Take B{ui,e) with 
e > and i = 1,2,3 three tiny balls and set /o and fi absolutely continuous 
with respect to the Lebesgue measure with densities lB{ui,e) + '^B(u2,e) and 
'^B{u2,e) + ^B{u3,e) respectively. The unique optimal transport plan for the 
usual quadratic cost is T2(x) = x + {u^ — Ui) if x E B{ui, e) and the identity 
elsewhere, but this map is not admissible for the cost c defined by (11.11) when 
C is the unit ball. The unique optimal transport plan for this cost is hence 
Tc{x) = X + (uj+i — Ui) if X G B{ui, e) for i = 1,2. 




Figure 1: The situation in Example 11.11 

As we said, the quadratic cost has often better properties than other 
general cost functions, and it is not astonishing that the first existence result 
for (M) was actually proven for c{x, y) = |y— a;p. Actually, Y. Brenier proved 
in 1991, existence and uniqueness of an optimal map T, and the crucial point 
in the proof is the existence of a dual formulation 



sup <^ / v{y) dfi{y) - / u{x) dfo{x) : v{y) ~ u{x) < c{x,y) 

where the sup is actually achieved. This allowed Brenier to use the primal- 
dual optimality condition and prove the optimality of the map 

T(x) = X — -Vm(x). 
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This scheme of proof was successfully extended to the cases c{x, y) = \x — y\^ 
with p > 1 (or, more generally, to any strictly convex function of x — by 
W. Gangbo and R.J McCann [13], and by L.A Caffarelli [6]. The original 
problem by Monge {p = 1) happens to be more difficult and needed other ar- 
guments and techniques. It was solved later (see L.C. Evans and W. Gangbo 
[TT] . N.S. Trudinger and X.J. Wang [16j, L. Ambrosio and A. Pratelli |3], M. 
Feldman R.J. McCann [12] )• The existence for the dual problem was also 
crucial in these proofs. 

When we add a constraint, the existence of an optimal pair of functions 
for the dual formulation [D) is not guaranteed: the usual proofs of this re- 
sult require continuity or finiteness of the cost (see the books of C. Villani 
|17j and [TS]). Some very irregular pairs of solutions (not even functions) 
were found for non finite Borel costs by M. Beiglbock, C. Leonard and W. 
Schachermayer Let us mention also that under the additional restrictive 
assumption J c{x,y) dfo{x)dfi{y) < +oo, W. Schachermayer and J. Teich- 
mann proved the existence of an optimal pair of Borel functions for (D) [T5] . 
Hence, though very interesting, these results cannot be used in our setting. 

Nevertheless, more recently, T. Champion, L. De Pascale and P. Juti- 
nenn developed a completely different idea and manage to solve the " L°° 
case" (see [ID])- It is the case where, instead of minimizing the integral 
J \T[x) — xl^dfolx), i.e. the norm, they considered the minimization of 
the L°° norm (as we mentioned before). This frame is very peculiar be- 
cause, due to its non-integral formulation, there is no dual formulation for 
the Monge problem. Hence, they developed a completely different idea, which 
was powerful enough to be extended later on to some cases where the cost is 
still integral, but for a function c(x, y) which is not strictly convex in x — y. 
In particular, this gave the proof of the existence of a solution T to Monge's 
Problem for c{x,y) = \\x — y\\ for a general norm ([8J and [9J), which was 
unsolved before. 

In this paper, we show that we can solve (M) with c defined by (11. ip using 
the same tools as T. Champion, L. De Pascale and P. Jutinenn. 

In the next section, we state the main results. In section 3, we make 
some comments and give a proof in dimension 1. The general proof is given 
in subsection 4.2. after proving some preliminary results in subsection 4.1. 
Finally, in the last section, we give an application of our result to the case of 
crystalline norms and to some L°° problems. 
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2 Main result 



In all the following, the cost c is defined by (11 .ip , /o and fi are two probability 
measures on R'^ and /o is absolutely continuous with respect to the Lebesgue 
measure. Let us introduce the relaxed formulation of the Monge problem 
introduced by Kantorovich. This new problem, (MK), satisfies inf(Mi^') = 
inf(M) (see [H] and [IB]): 

(MK) inf I / c(x, y) dn{x, y) : n e n(/o, A) 1 

where n(/o, /i) is the set of probability measure on M'^ x R'^ with fixed 
marginals /o and /i. The following classical result holds: 

Proposition 2.1. If there exists vr G Il{x,y) such that J^dy^-^d c{x,y) d7i{x,y) 
is finite, then (MK) has a solution. 

The classical strategy to solve (M) is to show that at least a solution 7 
of (MK) is concentrated on the graph of a map T i.e: 

J ip{x,y) d^{x,y) = J ip{x,T{x)) dfo{x) G CM')- 

Indeed, the map T is admissible for (M) if 7 is admissible for (MK) and 
it is optimal if 7 is optimal. Conversely for any T admissible for (M), 7 = 
(Id X T)jj/o is admissible for (MK), and if T is optimal, 7 is optimal too. 

In order to state the main theorem of this paper, we need to give some 
definitions on the geometry of convex sets. 

Definition 2.2. Consider a convex set C C M'^. 

1. If Xo is any point of C , we call dimension of C the dimension of 
Span(C — Xo), i.e. the dimension of the smallest affine space containing 
C . This number is denoted dim(C) and does not depend on the choice of 
Xq. The interior and boundary of C in the canonical topology associated 
to the Euclidian distance in such an affine space are called relative 
interior and relative boundary and denoted by liC and idC. 

2. A convex subset F G C containing more than one point is called a 
straight part of C if it is contained in relative boundary idC . 

3. A subset F of a convex set C is called a maximal fiat part of C if it 

is maximal for the inclusion among straight parts of C. 



5 



4- A subset F of a convex set C is called a flat part of C if it there is a 
finite chain of convex sets F = Fi C F2 C ■ ■ ■ C F^-i C F^ = C (k > 
2 ) such that, for each i = l,...,k — l,Fiisa maximal fiat part of Fj+i . 
Notice that in this case one has dim(Fi) < dim(F2) < ■ ■ ■ < dim{Fk). 
In particular, for every flat part F of C we have 1 < dim(F) < (d—l). 

5. A set C dW^ is said to he strictly convex if it has no flat parts. 

The aim of tliis paper is to prove the following result: 

Theorem 2.3. Assume that: 

i) there exists n G n(/o, /i) such that /j^d^j^d c(x, dn^x.y) < +00, 

ii) C has at most a countable number of flat parts. 

Then, if /o is absolutely continuous with respect to the Lebesgue measure on 
M.'^, any solution 7 of (MK) is concentrated on the graph of a map T which 
is optimal for (M) . 

Remark 2.4. Assumption ii) is obviously satisfied if C is strictly convex. 
Moreover, if C satisfies ii), then any flat part of C is a convex subset of an 
affine space M'^ for k < d which satisfies the same assumption. Finally, this 
assumption allows every bounded convex set in (because flat parts should 
be subsets of positive length of its boundary, which has finite perimeter, and 
hence they cannot be more than countably many) and all compact polytopes 
in any dimension. 

Corollary 2.5. Under the assumption of Theorem \2.3l the solutions of 
(MK) and (M) are unique. 

Proof of the corollary. The idea is classical : let 71 and 72 be two solutions 
of {MK), let Ti and T2 the associated transport maps (given by Theorem 
I2.3p . The probability measure 2i±22 jg qI^q associated to a transport map T3. 
Then for any G Cc(M'^ x R'^): 

1/2 ( / (p{x,Ti{x)) dfo{x) + ip{x,T2{x)) dfo{x)] = f{x,T3{x))dfo{x). 

This implies fo({x : Ti{x) 7^ T2{x)}) = 0, that is 71 = 72 and Ti = T2 almost 
everywhere. □ 
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3 Proof and comments in dimension 1 



The one-dimensional case is a very special one. Actually, for a very wide class 
of transport costs, the optimal transport is always characterized in the same 
way and is always the same, namely the unique nondecreasing transport map 
sending /q to /i. This is summarized in the following theorem (the proof is 
written here in its full generality, even if it is essentially taken from [T]): 

Theorem 3.1. Let : M — )■ R U {+00} be a l.s.c. strictly convex function 
(which actually means strictly convex on its domain {0 < +00}), and /q 
and fi two probability measures on M, with /o non-atomic. Then the optimal 
transport problem 

inf jy ^ 0(2/ - rf7r(x, : 7rGn(/o,/i)| 

has a unique solution, which is given by •y = (Id x T)#/o, where T is the 
nondecreasing transport map from /q to fi, which is well defined /q — a.e. 

Moreover, if the strict convexity assumption is dropped and is only con- 
vex, then the same T is actually an optimal transport map, but no uniqueness 
is guaranteed anymore. 

Proof. We will use the usual strategy based on c— cyclical monotonicity, i.e. 
the fact that any optimal 7 is concentrated on a cyclical monotone set F, 
where we can also assume (f){y — x) < -\-oo. This means that (x, y), (x', y') G F 
implies 

0(y - x) + (Piy' - x') < (Piy' -x) + (t){y - x'). (3.1) 

We only need to show that this implies (in the strictly convex case) a mono- 
tone behavior : we will actually deduce from fl3.ll) that x < x' implies y < y'. 
This allows to say that F is included in the graph of a monotone (nonde- 
creasing) multifunction. Hence, up to a countable set, for every x there is 
unique image y = T{x), and the map T is nondecreasing. 

To prove y < y' suppose by contradiction y > y' and denote a = y — x, b = 
y'—x' and c = x'—x. Condition (13.11) reads as 0(a)+0(6) < 0(6+c)+0(a— c). 
Moreover, we have 

6 + c = (1 — t)6 + ta, a — c = tb -\- (1 — t)a, for t = -. 

a — 

The assumption y' < y reads as 6 + c < a, which gives t g]0, 1[ (since it 
implies b < a, and hence t > 0, and c < a — b, and hence t < 1). Thus, 
convexity yields 

0(a) + 0(6) < 0(6 + c) + 0(a - c) 

< (1 - t)0(6) + t0(a) + t0(6) + (1 - t)0(a) = 0(a) + 0(6). 
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Since we assumed (pi^a) , (pip) < +00, we also deduce (pip + c),0(a — c) < 
+00 and the strict convexity implies a strict inequality, thus getting to a 
contradiction. 

This shows that a strictly convex cost always admits a unique optimal 
transport plan, which is given by the monotone nondecreasing transport map. 
The statement when is only convex is obtained by approximation through 

— x) + — (or through any other strictly convex approximation). □ 

It is important to notice that in this one-dimensional convex framework 
the optimal transport does not really depend on the cost. This is very useful 
for approximation procedures (if, for instance, one considers costs of the form 
\y — x\^ and lets p — )■ cxd it is possible to deduce the optimality of the same 
nondecreasing transport for the L°° transportation problem, see again jlOj). 
Not only, in our case when 0(z) = |zp for z G C and = +00 outside C, 
this shows that the optimal transport is not really affected by the constraint 
(differently from the example in that we gave in the introduction). What is 
affected by the constraint is the possibility of transporting /q onto /i at finite 
cost. Depending on the two measures, it is possible that the nondecreasing 
map T satisfies T{x) — x G C or not. If yes, then it is optimal; if not, then 
the problem has no solution with a finite cost. 



4 Proof of Theorem 2.3 



4.1 Preliminary results 

Let 7 be an optimal transport map for {MK). We recall that when the cost 
is l.s.c. (which is the case of our cost function c), then any optimal 7 is 
concentrated on a c-cyclically monotone set F. This means in particular 

(xo, 1/0), (xi, yi) eV ^ c(xo, Vo) + c(xi, yi) < c(xi, yo) + c(xo, yi) 

(we refer again to [T] for the actual definition of c— cyclically monotone sets, 
which imposes a similar condition for families of k points in F and arbitrary 
permutations of their second coordinates; here we only need the condition 
for k = 2). In our case, thanks to the definition of c, if we use the equivalence 

\xo - yo\'^ + |xi -yi\^ < |xo - yi\^ + \xi - yoP <^ {xi -xq, yi-yo) >0 
then we get 

If (xo, 1/0), (xi, 2/1) e F are such that: 1 , x\>n ^4 1^ 

2/o-xieCand?/i-xoeC > ^ {Vi ' Vo^Xi - Xo) > 0. (4.1) 
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Without loss of generality, we may assume that F is a-compact (which is 
useful for measurability reasons) and that any {x, G F is such that y — x E 
C. 

Arguing as in we introduce the following object: 

Definition 4.1. For any r > and any y G W^, we set: 

V-'(B{^)) := 7r,.(F n x 

where B{y,r) denotes the open ball of radius r > centered at y. The as- 
sumption of a— compactness on F guarantees that the set that we just defined 
is measurable. 

We recall the following results extracted from [10] and [8]: 

Lemma 4.2. The measure 7 is concentrated on a a-compact set R{T) such 
that for all {xq, yo) in R{T) and for all r > 0, the point xq is a Lebesgue point 
ofT-\B{yo,r)), t.e.: 

y \T~\B{y^,r))nB{xo,e)\ 

hm 1 — -] = 1 . 

£^0 \B{xQ,e)\ 

Proposition 4.3. Let (xq, Vq) G -R(F), let r > 0, 5 g]0, 1[ and ^ a unit vector 
in W'-, then for any e > the following set has a positive Lebesgue measure: 

v-\'B{^) n i?(xo, e) n C(xo, 5), 

where the set C(xo,^,5) is the following convex cone: 

{x: (x-XcO > {I - 5)\x - Xo\} 
(notice that it is an open cone, with Xq ^ C{xq,^,5)). 

Lemma 4.4. Let x E C , a, unit vector rj > such that x + t]^ E Int(C). 
Let < rj' < rj be fixed. Then there exist r > and 6 g]0, 1[ such that: 

\fyeB{x,r)nC, (c{y,^,6) HBi^) Clnt(C). 
The direction ^ is called an entering direction at the point x. 

Proof. First of all notice that, if x + r)^ E Int(C), then all the interior points 
of the segment connecting x to x + r]^ must belong to Int(C). 



9 



We prove the existence of r > and and 6 g]0, 1[ such that: 



yyeB{x,r)nC, (Ciy,^,6)n B{y,i)\B{y,r,'/2) C Int(C). (4.2 



If we prove (14.21) . the result then follows by passing to the convex envelope. 
The thesis (14. 2 p may be proven by contradiction. Suppose that it is false, 
and find a sequence ?/„ G C such that ?/„—)■ x, together with a sequence 

Zn i Int(C) and z„ G (c(?/„,^,5„) n B{yn,r]')\B{yn,i /2) for 5„ ^ 0. 

Up to subsequence we may suppose Zn — > z. This limit point z must satisfy 
z ^ Int(C) (since Int(C) is open) and |z — x| G [r]' /2,ri']. Moreover, from 
{zn-yn,0 > {l-Sn)\zn-yn\, we get (z-x,^) > \z - x\, i.e. 2;-a; = for 
A > 0. Then, we get A G [rj'/2,ri'] and finally we obtain a contradiction since 
z belongs to the interior of the segment connecting a; to a; + r]^, but not to 
Int(C). □ 



4.2 The proof in dimension d > 1 

We start proving a fundamental lemma: 

Lemma 4.5. Let (xo,?/o) <ind (xo,?/o) two points of R(T). We assume 

{(1 - t){yo - xo) + t{y', - xo) : t G [0, 1]} n Int(C) ^ 0. 
Then yo = y'o- 

Proof. To show this lemma, we proceed by contradiction and adapt the 
method developed in |10j to our case. Let {xo,yo) and (xo,?/o) two points 
of -R(r) satisfying the lemma's assumption with yo ^ y'^. 
We set ^ := j^^z^- The direction ^ is an entering direction at y^ — XQ and — ^ 
is an entering direction at y'^ — Xq. Using lemma U31 we take r > 0, p g]0, 1[ 
and e > 0, such that: 

(c(y[,-xo,-e,p)n5(y^-xo,£)) Clnt(C), (4.3) 

G 5(1/0 -Xo,r + £), (c{z,i,p)nW^^ Clnt(C). (4.4) 
Moreover, if p and r are small enough, we also have 

Vx G C(xo,e,p), Vy G 5(yo,r) : (x - Xq, y - Z/q) < 0. (4.5) 

Indeed the direction of x — xq is almost that of (up to an error p) and that of 
y — y'o almost — ^ (up to an error of order r/|?/o — 2/ol)- The two vectors being 
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opposite, it is clear that if both 6 and r are small then the scalar product 
stays negative. By use of Proposition 14.3^ we get the existence of a couple 
(x, G r such that x E B{xo, e) fl C{xo, C,, p) and y G -B(|/o, f)- 
Then the couple (x, y) eT satisfies: 

• y'^ — X = y'^ — Xq + {xq — x) G C, as a consequence of fl4.3p . 

• (x-Xo,i/-i/o) < 0, by (gSl). 

To get a contradiction with (14 .ID . it remains only to show y — xq E C. This 
comes from (14.41) . Indeed y — x = {yo — xq) + {xq — x) + {y — yo) is in 
B{yo — Xo,r + e) and x — Xq E (C(0, ^, p) fl B{0, e)) so that 

y — Xq = {y — x) + {x — Xo) E IntC. □ 

We easily deduce 

Corollary 4.6. // C is strictly convex, then any solution of (MK) is sup- 
ported on a graph. 

Moreover, if C is not strictly convex, any solution of {MK) is supported on 
a set R(T) satisfying the following property: 

If {xoyUo) (^nd {xojy'o) are in R(T), then (yo — Xq) and (?/q — Xq) belong to a 
same flat part of C. 

To prove Theorem 12.31 itself, we proceed by induction on the dimension d 
of the space. The theorem in dimension d = 1 has already been proved. Let 
us show that if the theorem holds in any dimension lower or equal to d — 1, 
it also holds in dimension d. 

Denote by Co the set of points of C which do not belong to any flat part of 
C, and by (Cj)jg/ the family of the maximal fiat parts of C. Let 7 an optimal 
transport plan and R(T) defined as previously so that 7 is concentrated on 
RiT). We set 

Ti:= {{x,y) : y-xE d} H RiT), G / U {0}. 

By Lemma 14. 5^ the measure 7 L Fq is concentrated on a graph, and, for any 
i E I, we have the following property: 

(x, y)ET,^ ({(x, y') : (x, y') E R{T)} C T,). 

Then it is enough to show that for any i E I, 'ji := 7Lrj is concentrated on a 
graph and to glue the results together. Each measure 7^ is an optimal trans- 
port plan for the same cost c between its marginals. The assumption that 
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/ is at most countable allows to decompose 7 into a sum of these measures, 
thus avoiding the need for disintegration results. In particular, it is clear that 
the key assumption that /o is absolutely continuous stays true for the first 
marginal of 7^ (which would not be easy in case of a more-than-countable 
disintegration) . 

Let us fix i e 7, let n < d — 1 be the dimension of Cj and z a fixed vector in 
Cj. Up to a change of variables, we may assume 

Span(Ci — z) = Span(e^, ....e"") 

where (e^, .... e^) is the canonic base of M^. 
Notations: We denote: 

7 := 7Lri, /o := -K^l^, fi := TTy^. 
Note that y — x E Ci imphes y2 — ^2 — Z2 so 7 is concentrated on 

{{x, (yi, X2 + Z2)) : x = (xi, X2) G yi G W"}. 

Making a disintegration of 7, /o, /i we write: 

7(0;, (yi, X2 + Z2)) = fi^'^i^i, yi) ® C(x2, y2), with C := Tr^2,y2^1, 

/o = /^^(xi) ® (7r.,tt/o) (X2), A = /r(yi) ® (nyM) (y2). 

Notice moreover that the measure ( is concentrated on the pairs of the form 
(x2, X2 + Z2). Then we have: 

J c{x,y) dj{x,y) = y |xi - yil^ ci7(x, y) + |z2|^7(ri) 

|xi - yi|^ rf/i''2(xi,yi)dC(x2,X2 + Z2) + |z2|^7(rj). 

Weset: c(xi,yi) = | '^^-^^1' 'Mji -^1:^2) e Q, 
\ ujij \ +00 otherwise. 

Lemma 4.7. (- almost every (x2,X2 + Z2), we have: 

fo^ — ^''xitl/t''^, and ff^ = TTy^Jj/t^^, moreover /q^ is absolutely continu- 
ous with respect to the Lebesgue measure, 
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ii) the measure fi^'^ on M.'^ x Mf^ is a solution of the following transport 
problem: 

(P-^) minjy" c(xi,yi) d7r(xi,yi) : tt G n(/o''% . 

Then using the theorem 12.31 in dimension n < d — 1, we get that, for ^- 
almost every (x2,X2 + Z2) the measure fi^^ is concentrated on a graph. Then, 
recalhng that ( is concentrated on {(x2,X2 + Z2)}, we get the desired resuh. 

Proof of the lemma, i) Let us show /^^ = n^j^^fi^'^. On the one hand, we 
have: 

/o = 7rxi,x2tl7 = 7rxi,x2tl(/i'''(xi,yi) ® C(x2,y2)) 

= (7rx,tl/i''^)®(7rx,tlC(x2,y2)). 

On the other hand: 

7^X2 tic (X2,y2) = 7rx2tl7rx2,y2tl7(x,y) = TTxJ/o. 

Then, by uniqueness of the disintegration /^^ of /o outside a TTxatl/o neghgible 
set, we get the result. 

In the same way, we can show ff^ = vTyj^/i^^. Moreover as /o is abso- 
lutely continuous with respect to Lebesgue, also is /^^ (with respect to the 
ra— dimensional Lebesgue measure, instead of the rf— dimensional one), 
ii) We proceed by contradiction. Assume there exists E C {(x2,X2 + Z2)} not 
negligible such that, for all (x2, X2 + Z2) G E, the measure fi^^ is not optimal 
for {P^'^)- For any (x2,X2 + Z2) G E, let us pick another measure ft^^ chosen 
so as to be optimal for (P^^)- We set: 

^ = fi''^ (g) C\_E + fx""' ^ C^-E" 

where E'^ is the complementary of E in Ci. 
Then 7 G n(/o, /i), indeed: 

VTxtlT = vr4(/i''2®CL£;) + 7r4(/i"2®CL£;'=) 

= TTx, m^'' ® VTx, HC L ^ + VTx, H/i"^ ® TTx, K L E' 

= f^'^in^.C^E + n^JC^E^) 
= /o^®vrx,C = /o, 
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the same holds for the y marginal. Moreover: 

J c{x,y) d'y{x,y) = j j \xi - yi\'^ dp.'^^{xi)dC{x2,^2 + Z2) 

+ j j |xi-yip rf/i"^(xi)rfC(x2,X2 + Z2) + |z2p7(A 
< / / |xi - yil^ (i/i''2(xi)(iC(x2,X2 + Z2) 



'E 



xi - yil^ ci/i''2(xi)ciC(x2,X2 + Z2) + |z2p7(Ci 
c{x,y) d^{x,y). 

This is impossible because 7 is optimal for the transport problem with cost 
c between its marginals /o and /i. □ 

5 An application to crystalline norms and 
problems 

We present in this section how the results of this paper may be applied, 
following the strategy presented in to the problem of the existence of 
optimal transport maps when the cost is a crystalline norm ||x — y\\. We 
recall that a norm is called crystalline if its unit ball is a convex polytope 
containing in its interior with a finite number of faces; In particular, this 
means that it has the form 1 | = maxj z ■ Vi for a certain finite set of vectors 
(t>j)j^]^ jj. This case has been first studied in [2] and was one of the first 
steps towards the extension of the Monge problem to general norms. The 
existence of an optimal map for norm costs has finally been established in 
full generality in 

The strategy proposed in [7J concerns the minimization of a transport 
cost j c{x — y)d7T for a convex, but not necessarily strictly convex, function 
c. In the same spirit of the proof of Theorem 12.31 a decomposition may be 
performed on an optimal transport plan 7, according to the "faces" of the 
cost. In this case, since the crystalline norm has a finite number of faces 
Fi (which are the cones on the maximal fiat parts of its unit ball), one can 
write 7 = 7j where the measures 7* are obtained in the following way. If 
fo« CJ^ ■, duality implies (see [7]) that for there is an index i = i(x) 

such that {x,y) G spt(7) implies x — yEFi{m case of non- uniqueness of this 
index just pick one at random, for instance by fixing an order relation on the 
faces). Then define 7i = 'y\-{{x,y) : i{x) = i}. 
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This means that it is enough to prove that every plan 7* may be actually 
induced by a transport map (more precisely, that there exists a new plan, 
with the same marginals and an improved transport cost, which is induced 
by a map). This is quite easy since the cost, when restricted to a face, is 
actually affine, and hence every transport plan with the same marginals and 
supported on {{x,y) : x — y E Fj} gives the same cost. The issue is thus 
reduced to finding an arbitrary admissible transport map T, with prescribed 
marginals /g and /{ (defined as the marginals of the original 7*), satisfying 
X - T{x) e Fi. 

Thanks to the results of the present paper this is possible by considering, 
for instance, the minimization problem (MK) with convex set C = F^. This 
selects a particular transport map (the one minimizing the quadratic cost 
among admissible ones), but gives in particular the existence of one such 
transport map. The only point to be noticed is the fact that a transport 
plan with finite cost is actually known to exist, and it is 7*. 

As a consequence, we have a new proof of 

Theorem 5.1. If fo << dx and || ■ || is a crystalline norm, the transport 
problem 

minjy \\x-y\\dTX : 7rGn(/o,/i)| 

admits at least a solution induced by a transport map T . 

Analogously, the arguments used in this paper may also be applied to 
problems like those studied in [TD]. We already pointed out that the strategy 
of the quadratic perturbation is useful when one only wants to prove the 
existence of at least a transport map T sending /o onto /i and satisfying 
T[x) — X E C. Suppose now that we want to solve 

(Moo) min{7r - esssup II?/ - x|| : vr G n(/o, /i)} , 

where 1 1 ■ 1 1 is a (possibly asymmetric) norm whose unit ball is C. Notice that 
we have the equality vr — ess sup I ly — x| I = max{| |y — a;| I : {x,y) E spt{n)} 
and that this quantity measures the maximal displacement in terms of the 
norm 1 1 ■ 1 1 . Then it is sufficient to denote by L the minimimal value of such 
an ess sup and to notice that tt is optimal if and only if it is concentrated on 
{{x,y) : y — X E LC}. In order to find a solution induced by a transport 
map T it is enough, for instance, to take a map T solving (M) for the cost c 
defined with the dilated body LC instead of C. Hence we have the following 
theorem, generalizing the results in (and using, by the way, a different 
selection principle, since the plan which was selected in [10] was the infinitely 
cyclically monotone one, and not that minimizing the quadratic cost). 
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Theorem 5.2. Suppose that C satisfies the assumptions of Theorem \2.3[ 
that G int(C) and that /g << dx. Let \ \ ■ \ \ be defined through the formula 
\\z\\ := inf{A > : z/X G C}. Then, the L°° transport problem (Moo) admits 
at least a solution induced by a transport map T. 
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